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We consider the pattern dynamics of the lamellar phases observed in Rayleigh-Benard convec- 
tion, as described by the Swift -Hohenberg equation, and in the weak segregation regime of diblock 
copolymers. Both numerical and analytical investigations show that the dynamical growth of the 
characteristic length scale in both systems is described by the same growth exponents, thus sug- 
gesting that both systems are members of the same universality class. 
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I. INTRODUCTION 

The study of the dynamics of pattern formation in sys- 
tems far from equilibrium encompasses examples from 
both physics, chemistry and biology Despite com- 
pletely different physical origins some systems exhibit 
identical morphologies and pattern dynamics, and may 
be perceived as members of the same universality class. 

In this paper we consider the pattern dynamics of 
two morphologically identical systems, namely Rayleigh- 
Benard convective rolls and weakly segregated diblock 
copolymers. At short times after a quench from the uni- 
form stable phase to the unstable phase both systems 
develop a labyrinthine domain morphology consisting of 
rolls (or lamellae) of a well-defined width w. Initially the 
rolls are randomly oriented, but as time increases they 
locally align up in parallel thereby creating an increas- 
ingly ordered pattern [Fig.El. We have investigated the 
dynamics of this coarsening process by numerical integra- 
tion of the appropriate Langevin equations and by ana- 
lytical considerations. Both approaches agree that the 
characteristic length scale of the systems scales dynami- 
cally with growth exponents which are common to both 
systems, thereby suggesting that the pattern dynamics 
of Rayleigh-Benard convection and diblock copolymers 
belong to the same universality class. 

The observed ordering phenomenon is driven by two 
mechanisms, namely interface relaxation and defect an- 
nihilation. The effect of the former mechanism can in a 
defect-free system be calculated by considering the speed 
at which a modulated interface relaxes to its (straight) 
ground state. Specifically we apply the projection op- 
erator method Q developed for interface relaxation in 
the Rayleigh-Benard system to the same problem in di- 
block copolymers, thus providing a systematic treatment 
of both systems. Furthermore we show how the applica- 
tion of a general approach to interface relaxation recently 
developed by one of the authors H leads to the same re- 
sult for the Rayleigh-Benard system. 

This paper is organized as follows: In Sec|n] we intro- 



duce the two models we study. Our numerical work is 
presented and discussed in Sec.|l|. Secjl^ contains the 
theoretical considerations including a brief review of the 
projection operator method, and Secj^ concludes with a 
summary and discussion. 




FIG. 1. The coarsening process. The figure shows snap- 
shots of the domain configurations in the diblock polymer 
system shortly after the quench from the disordered to the 
bistable phase (a) and at increasingly later times (b)-(d). The 
pictures are contour plots of 128 x 128 systems where the con- 
tours are defined by 4>(x, y, t) = 0. The order parameter field, 
cf>(x,y,t), was obtained by numerical integration of Eq.(^) at 
zero thermal noise. Simulations of the Swift -Hohenberg sys- 
tem [Eq.(^)] produces domain configurations which morpho- 
logically are indistinguishable from those here presented for 
the diblock copolymer system. 



II. MODELS 

In the Rayleigh-Benard system a simple fluid is con- 
fined between two horizontal plates which are heated 
from below, and for values of the Rayleigh number, R, 
larger than a critical value, R c , an instability occurs 



1 



which transform the uniform state to a state consisting 
of spatially periodic convective rolls. Near the onset of 
the convective instability the free energy functional, F, 
of the Raylcigh Bcnard system is, in dimcnsionlcss vari- 
ables, well approximated by the form 



F\ 



dV -(fc o 2 +V 2 ) 2 ]^/2 + 4 /4), (1) 



deduced by Swift and Hohenberg (SH) [Q. Here the 
scalar order-parameter field <f> — <f>(x,y,t) is related to 
the local vertical fluid velocity, e = (R — R c )/R c is the 
reduced Rayleigh number which acts as the control pa- 
rameter of the system, and fco is the wavenumber corre- 
sponding to the period, A = 2w, of the modulated struc- 
ture, i.e., fco = tt/w. 

For small e the order parameter field is locally well 
described by a single mode approximation, <^>(r) ~ cos(k- 
r), where k is perpendicular to the orientation of the rolls, 
and as e — * this approximation is exact || ■ Minimizing 
the free energy Eq.(Q) in the single mode approximation 
yields fc = ko as the selected wavenumber of the steady 
state. As customary we use ko = 1. The equation of 
motion for <p is given by the Langevin equation dt<j> = 
—SF[4>]/6<p + C, where F is the above free energy and £ = 
C(r,t) is thermal noise correlated as (£(r, t)£(r', t')) = 
2A5(r — r')5(t — t'), where A parameterizes the strength 
of the thermal fluctuations. Thus the Swift-Hohenberg 
equation reads 



d t ^ = e(t>-{kl + v 2 ) 2 4>-4? + C- 



(2) 



A diblock copolymer (DC) is a linear chain molecule 
joined together by two strings of equal length of e.g. A 
and B monomers. The polymerization index, N, is thus 
N = Na + Nb where Na — Nb are the numbers of A 
and B monomers, respectively. Above the critical tem- 
perature T c , A and B mix, whereas below T c the two se- 
quences are incompatible and the copolymer melt under- 
goes phase separation. However, spinodal decomposition 
U cannot continue indefinitely because of the chemical 
bond between the sequences. As a result the phase sep- 
aration occurs on a length scale bounded above by the 
length of a stretched polymer chain (typically less than 1 
micrometer) where banded domains of A-rich and £>-rich 
regions alternate in the final equilibrium state. The free 
energy of a diblock copolymer melt below T c is given (also 
in dimensionless variables) by a modified Cahn-Hilliard 
free energy functional 



F[<p] 



d d r [f(4>) + (l/2)(V</>) 2 



(r/2) / d d rd d r^(r)G(r,r>(r') 



(3) 



where <fi(r,t) = </>a(i", t) — 0s(r,t) is the local con- 
centration difference between the A and B monomers, 
f((f>) = — </> 2 /2 + </> 4 /4, is the bulk free energy density, 
and r is a control parameter inversely proportional to the 



square of the polymerization index, F ~ 1/N 2 . Finally, 
the Green's function, G(r, r'), in the second integral is de- 
fined by the Poisson equation, V 2 G(r,r') = — S(r — r'). 
The order parameter for this system is a conserved quan- 
tity, thus the appropriate Langevin equation for the time 
evolution of </> subsequent to a quench from the disor- 
dered to the bistable phase, is dt(j> — V 2 8F[<fi]/S<fr + C, or 
inserting Eq.(||), 



V 2 (- 



(4) 



where the noise £, representing the effect of ther- 
mal fluctuations, has the correlations (C(r, t)((r' , t')) = 
-2A\7 2 5(r - r')S(t - t'). For T just below the critical 
value, r c = 1/4, Eq.(^) describes the dynamics of weakly 
segregated lamellar domains with a well-defined width 
w = n/ko, where fco = T 1 / 4 is the wavenumber which 
minimizes the free energy [Eq. ([}])] in a single mode ap- 
proximation ||. 

The diblock copolymer equation [Eq.Q] can conve- 
niently be rewritten in a form resembling the Swift- 
Hohenberg equation [Eq.(S)], 



dt<j) = ecf>- (1/2 + VT + V>' 5 + C 



(5) 



where e = T c — T. Linearizing in Fourier space about 
(f> = we find, in both Eq.(||) and Eq.(||), that fluctua- 
tions, 6<j>k, in the order parameter decay exponentially, 
S4>k(t) = S(j>k(0) exp[— uikt] with rate cut — (a — fc 2 ) 2 — e, 
where a = 1,1/2 for the SH and DC system, respec- 
tively. Thus both systems have a band of wavevectors, 
fc_ < fc < fc + , k± = ya±\/e, for which the uniform 
state is unstable. In the nomenclature of Cross and Ho- 
henberg [Q this means that both systems are stationary 
periodic, or Type I s . 



III. SIMULATIONS 

We have solved the SH and DC equations numerically 
using a finite difference scheme on two dimensional lat- 
tices of size 512 x 512, with periodic boundary conditions. 
Numerical algorithms for the spatio-temporal evolution 
of both systems were obtained by replacing, in Eq.(j^) 
and Eq.(|, dt<f>(r,t) by (<^ +1 - ^)/At, and V 2 0(r,i) 
by the following discretized Laplacian 



(Ax) 2 



2^1^ 10 

(NN) (NNN) 



"*3 > 



(6) 



which includes contributions from both nearest neigh- 
bours (NN) and next-nearest neighbours (NNN). Here 
the indices represent the coordinates (x, y) and the 
index n represents time. Connection to absolute time 
and spatial coordinates is established through the re- 
lationships t — nAt and r = (ix + jy)Ax. The spe- 
cific choice of coefficients in Eq.(ff) ensures that the 
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Laplacian, in Fourier space, is isotropic to second or- 
der in fc 2 , i.e., the form of the Fourier transform, 
r k , of Eq.(§), is T k = -fc 2 + const, x (Ax) 2 k 4 + 
0((Ax fc 2 ) 3 ). For the diblock copolymer system the 
fluctuation-dissipation relation for the discrete equation 
can be maintained by generating two independent Gaus- 
sian variables v;] (n), {n) with zero mean and corre- 



lations {v\f {m)v$ (n)) = 2AAt Si. k Sj,iS m , n S atb and then 
i r„(i) _„« , „(2) _„(2), 



letting fe(n) = -^M+i, ~ </ + " In 

the simpler case of the SH-equation, Cij( n ) is a Gaus- 
sian distributed field with zero mean and correlations 
(dj( m )Cki(n)} = 2AAt5 i>k Sj t iS mtn . 

An inherent complication in this type of numerical sim- 
ulations is the conflicting constraints which the choice of 
the stepsizes is subject to. The need for numerical accu- 
racy requires (Ax, At) to be vanishingly small, whereas 
the finite computational power available requires the op- 
posite. Specifically, a linear stability analysis Jn| of 
the above algorithm with the Laplacian given by Eq.(j|) 
shows that, in order to avoid spurious solutions aris- 
ing from the subharmonic bifurcation, the dimcnsionlcss 
mesh size Ax and timestep At must satisfy the rela- 
tion At < 2/[(q - 16/[3(Ax) 2 ]) 2 - e], where, as before, 
a = 1, 1/2 for the SH and DC systems, respectively. In 
practice, the size of Ax is dictated by the the smallest 
length scale in the problem, which is the selected wave- 
length A = 2w. In order to avoid lattice pinning it is de- 
sirable to have many lattice points per wavelength. This 
quantity is given by A/ Ax, so by lowering Ax any number 
can be obtained. However, from the above stability re- 
lation we see that decreasing Ax below unity drastically 
reduces the maximum allowable size of the timestep, and 
hence increases the required computer time. 

We have performed our simulations using the values 
(e, Ax, At) = (0.25, 2tt/8, 0.025) for the Swift-Hohenberg 
system and (e, Ax, At) = (0.05, 1.0, 0.05) for the diblock 
copolymer system, where both set of values satisfy the 
appropriate stability relations. In the SH system the se- 
lected wavelength is approximately 2ir, so Ax = 2n/8 
gives 8 lattice points per wavelength. The corresponding 
quantity in the DC system, which we consider in the weak 
segregation limit or small e, is approximately 9, since here 
the selected wavelength is A = 2tt / (1/4 — e) 1 / 4 . 

Appropriate to a critical quench the systems were ini- 
tially prepared in the homogeneous single phase state by 
assigning to each lattice site a small random number uni- 
formly distributed about = 0. Nonzero temperatures 
were simulated using the fluctuation strengths A = 0.4 
and A = 0.1 for the SH and DC systems, respectively. 
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FIG. 2. Test of the scaling form Eq.(Q) and (inserted) time 
evolution of the structure factor illustrated with data from 
simulations of the diblock copolymer system at zero thermal 
noise (Here depicted in arbitrary units). The scaling collapse 
was obtained with the value x = 1/5 of the scaling expo- 
nent. The data sets {<>,*,+} represent the (dimensionless) 
times {1.8 x 10 4 ,5.6 x 10 4 ,1.8 x 10 5 }. 



A. Dynamical scaling 



We monitor the coarsening phenomenon by means of 
the usual structure factor S(k,t) — |</>k(i)| 2 , where </>k(i) 
is the Fourier transform of the order parameter. The 
circularly averaged structure factor, 5(fc,t), is sharply 
peaked around the wavevector ko which corresponds to 
the width of the rolls, and as time evolves it becomes 
increasingly sharper and higher. Assuming dynamical 
scaling, the simplest scaling form for the structure factor 
is 



S(k,t)=t x f(t x [k-k ]), 



(7) 



where f(y) is a scaling function. This form implies that 
the width, Afc, of the structure factor and its intensity, 
S(ko,t), scale as Afc ~ t~ x and S(ko,t) ~ t x . In agree- 
ment with previous work by a number of authors, our 
data from the SH system satisfies this scaling form with 
the scaling exponents x = 1/5 and x = 1/4 at zero and 
non-zero thermal noise respectively 1 1 . 1^] . Further- 
more, we find that the diblock copolymer system also 
obeys Eq.(Q) with the same values of the scaling expo- 
nents [Fig.(0)]. 
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FIG. 3. The local director field, n(r) = V0(r)/|V0(r)|, 
here illustrated as small bars, from which the correlation func- 
tion Eq.(^) is computed. For visual clarity only directors near 
the domain boundaries (solid contours) are depicted. 

A more direct method of probing the rolls increasingly 
orientational order is computing a correlation function, 
C nn (r, t), of the 'nematic' order parameter n = V0/|V^|, 
i.e. the unit vector normal to surfaces of constant <f> 
[Fig.(||)]. Explicitly we have computed the correlation 
function 

C nn (r, t) = 2- £([n(x + r, t) ■ n(x, <)] 2 ) - 1, (8) 

X 

where N 2 is the volume of the system and (. . .) means 
a statistical average implemented through several inde- 
pendent runs. We compute ([n(ri) • n(r2)] 2 ) rather than 
(n(ri) ■ n(r2)) since we are interested only in the rela- 
tive angle, 6>(n, r 2 ), between the directors at sites i"i and 
r 2 . For sites separated by large distances the correspond- 
ing directors can be expected to completely decorrelated, 
thus 2 (cos 2 e) e -l=0. 

The time complexity of the algorithm for both C nn 
and its corresponding structure factor, S nn Qs.,t) = 
N- 2 J2r Cn«(r, t) exp(ik-r), is N 4 , and with N = 512 ex- 
cessive computer time is demanded. This problem can be 
circumvented by introducing the two-dimensional tensor 
Qab(r,t) = n a (r,t)n b (r,t), where n a , a = {x,y} are the 
components of n. In terms of Q a b Eq.(||) appears as 

C m (r, t) = j^ 5I(g ab (x + r, t)Q o6 (x, t)) - 1, (9) 

X 

where summation over repeated indices is under- 
stood. Since S nn now has the form S nn (k.,t) = 
2(QabQL,t)Qab(—Kt)) ~ $k,o, C nn can quickly be com- 
puted via S nn using Fast Fourier Transform J13|. 
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FIG. 4. Time evolution of the director correlation function 
Eq.(^) illustrated with four successive (dimensionless) times 
{5.6 x 10 3 ,1.8 x 10 4 ,5.6 x 10 4 ,l-8 x 10 5 } increasing from left 
to right. We extract the time evolution of the length scale, 
L(t), by monitoring the r a (t) for which C n n{r a (t)) = a where 
a < 1 is some constant (the horizontal dotted lines show 
a = {0.2, 0.3, 0.4}). The scaling exponent y is extracted from 
a log- log plot of r a (t) versus t (inserted). The data shown 
result from a simulation of the SH system at zero thermal 
noise, and y assumes the value y = 0.25 ± 0.02. 

The expected scaling form for the correlation function 

C nn (r,t)=F(r/L(t)), (10) 

where F is a scaling function, and L is a length scaling 
as L(t) ~ t v . For both the Swift-Hohenberg and diblock 
copolymer system we find this scaling form to be satisfied 
with the scaling exponents y = 0.25 and y = 0.30 at zero 
and finite noise, respectively [Fig. ([!])]. 

The values of the scaling exponent y agree with the 
findings of Hou et al |l^] . These authors measure the den- 
sity, p(t) , of topological defects in the Swift-Hohenberg 
system and find the algebraic decay p(t) ~ t~ v where 
y = 0.25 and y = 0.30 at zero and finite noise, respec- 
tively. The boundaries between plane- wave domains con- 
sist of topological defects. Therefore the defect density 
must scale as the perimeter density of the domains which 
again scales as the reciprocal, of the linear size of 

the domains. Furthermore, Hou et al find that the en- 
ergy of the Swift-Hohenberg system, Eq. ([!]), decays as 
the defect density. Also here the diblock copolymer sys- 
tem behaves as the Swift-Hohenberg system. Measuring 
the energy, as given by Eq.(||), we find the algebraic de- 
cay E(t) ~ t~ v with the same values for y as above. 

IV. THEORY 

Theoretical analysis of the pattern dynamics of lamel- 
lar phases is complicated by the presence of topological 
defects and current theories apply only to systems with- 
out defects. However, locally Type I s systems exhibit 
nearly ideal lamellar structures where, in two dimensions, 
the order parameter can be described as an amplitude 
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modulated plane wave (p(r,t) — {<f)oA(x , y , t)e tk " x + c.c], 
where we have assumed lamellae perpendicular to the in- 
direction; A is a complex amplitude and c.c. denotes the 
complex conjugate. Inserting this form into the equations 
of motion, Eq.(Q) or Eq.(|^), we obtain in the absence of 
noise the Amplitude Equation 

T d t A = eA + e [d x - (i/2k )d 2 ] 2 A - g \A\ 2 A, (11) 

where To,£o and go are constants. The derivation of 
Eq.(|TTl) from the Swift-Hohenberg equation [Eq.(§)] is 
described in Ref. and the method of this reference 
can easily be extended to the DC-equation [Eq.(||)]. The 
amplitude equation describes the dynamics of both the 
magnitude, \A\, and the phase, 9(r,t), of the complex 
amplitude A. By perturbing the steady state solution of 
Eq.(O) we obtain, to lowest order in e, the Phase Equa- 
tion^ 



D\\d 2 x t 



D 



dp, 



(12) 



where Dm and D± are diffusion coefficients in the parallel 
and normal directions, respectively. Dimensional analy- 
sis of Eq.(|l2]) implies a i^-growth of the characteristic 
length scale, in disagreement with numerical investiga- 
tions which favour a smaller value of the growth expo- 
nent. However, as discussed below, by considering how a 
curved interface relaxes, working to second order in e, a 
transient regime with i^-growth can be predicted. 



A. Projection operator method 

In order to follow the slowly varying orientation of the 
rolls (or lamellae) Elder et al Q introduces a coordinate 
system that tracks the interface given by the points at 
which <j) = 0. Specifically the Cartesian coordinates x and 
y are mapped onto curvilinear coordinates (it, s), where u 
and s are locally normal and parallel to the lines </>(r, t) — 
0. Assuming that the curvature of the individual rolls are 
small, the Laplacian in the new coordinates becomes 



6? 
du 



d_ 

du 



ds 2 ' 



(13) 



where re is the local curvature. Assuming that the sta- 
tionary solution of the one-dimensional Swift-Hohenberg 
equation is a good approximation in the normal direction, 
Eq.(0) becomes 



d(j) s du 
du dt 



2k 



df_ 

du 



+ 



d 3 <p s 
du 3 



di 



(14) 



where k ss — d 2 K,jds 2 and 4> s is the solution of 
(j) s (u(r,t)) 3 = [e-(k 2 + dl) 2 ]cj) s (u(r,t)). The final term, 
A, in Eq.(p^) contains terms of higher order in re, and 
terms involving the derivative of re in the direction normal 
to the lamellae: A = (n uu + nn u )d u (j) s + (2re„ + K 2 )d 2 4> s , 
where k u means d u n, etc. Application of the projection 
operator, 



2tt 



du d u (, 



to Eq.([f4j) produces the final result 



(15) 



(16) 



where v = dtu 
-2(k 2 + (3/a), a = 



is the interface velocity, a 
(ko^J^lduidur) 2 , and 



(k /2n) J^° ko du{d u (j) s ){dl<j) s ). The term involving A 
drops out from the final result, because Ad u <fi s can be 
written as ^[(re^ + re 2 /2)(d u (j) s ) 2 ], which vanishes when 
integrated over one lamellar thickness. 

In order to evaluate the coefficient a in Eq. (|l^) the sta- 
tionary solution is expanded to leading order in e, yield- 
mg |. 



4> s (u) = $i cos(fc it) + $3 cos(3k u), 



(17) 



with coefficients $i = and $ 3 = — <&f /256. Using 

this expansion we find a = e 2 /256, remembering that 
k 2 = 1. 

Applying the same analysis as above to the DC- 
equation, the equation corresponding to Eq.(n4) becomes 



d(j> 3 du 
du dt 



df_ 

du 



.d 3 



du 3 



3(<n s 



du 



-Kss^— +A 

du 



(18) 



where s is the solution of d 2 4> s (u(r, t)) 3 = [e — (1/2 + 
d 2 ) 2 ](jf (u(r , t)), and A has the same meaning as in 
Eq. (|l4[) . Using the projection operator |Eq.(|l5|)l we re- 
trieve Eq. (|l~6|) , only now with a = — (1 + 2/3 /a — 

where 7 = (fc /27r) jl^° kQ du (cf) s ) 2 (d u (f> s ) 2 , and a and 
(3 are as defined above. The quantities u, (3, and 7 
can be determined by substituting the form Eq. ( |l7j ) into 
the free-energy functional Eq.(||), and minimizing with 
respect to fcoi *J?i , and $3. To the required order in 
e=l/4- T, the result is k = T 1 / 4 = (1/4 - e) 1 / 4 , 
$f = (8/3)e+ (19/6)e 2 , and $ 3 = -(9/128)$?, leading 
(after some algebra) to a = (45/32)e 2 , correct to leading 
non-trivial order in e. 

In this approximation a is a very small number, a ~ 
2.4 x 1CT 4 and a ~ 3.5 x 1(T 3 in the SH and DC systems 
when e = 0.25 and e = 0.05, respectively. Dimensional 
analysis of Eq.(fl6l) therefore implies a crossover in the 
growth of the characteristic length scale from a transient 
i 1 / 4 -growth to an asymptotic t 1 / 2 -growth. The crossover 
occurs approximately when (at) 1 / 2 — i 1 / 4 , that is, when 
t ~ 1.7 x 10 7 in the SH system and when t ~ 8 x 10 4 
in the DC system. These crossover times far exceed the 
latest times we have been able to probe in our simula- 
tions, but though (as pointed out by Elder et al an 
appealing interpretation of the numerical results is that 
they witness the transient regime, there is no numerical 
evidence of any crossover behaviour. However, since the 



5 



estimated crossover time for the DC system is three or- 
ders of magnitude smaller than the corresponding time in 

the SH system, the DC system is the obvious candidate (a) 

for future investigations. ^J k <> 



B. Relaxation of a modulated interface 

The same problem can be investigated using a general 
approach to growth exponents recently developed by one 
of the authors || . 

We consider a small regular perturbation of the per- 
fect lamellar phase and wish to determine the rate at 
which the system relaxes to its ground state. Setting 
<fi(x,y,t) = 4>L(x) + (f>(x,y,t), where 4>l{x) is the station- 
ary lamellar solution of the Swift -Hohenberg equation 
[Eq.<§], 



= e0 L - (fcg + dlf4> L - 4> 



3 

Li 



(19) 



and is a small perturbation, the linearized equation of 
motion for becomes 



S t = e0-(fc o 2 +V 2 ) 2 0- 30|0. 



(20) 



A modulation of the lamellar phase with wave vector 
q < fc is <f>(x,y,t) = §l\x + A(t)cos(qy)} ~ 4>l{x) + 
<p' L (x)A(t) cos(gy), where 4>' L means d x <f>L and the ampli- 
tude, A, of the modulation is assumed small compared 
to the lamellar spacing [Fig. (^j)] . More generally, we can 
write the modulated phase as 



(x, y, t) = cj> L (x) + A e cos(qy)ip q (x) 



(21) 



where the function ip q (x) has the same periodicity as 
<t>L, {4>l ~ cos(fcoa;),fco = 1), and satisfies the condi- 
tion ip q= a(x) — (f)' L (x), since q = represents a uniform 
translation of the lamellae. We have written the time de- 
pendence of the amplitude as A(t) = A e~ UJqt , describing 
the decay of the eigen-perturbation of wave vector q at 
a characteristic rate cu q . 

Insertion of Eq.(pl|) into Eq.(|2(]) yields an eigenvalue 
equation, Hipq = ^qi^qi f°r the relaxation rate with 
'Hamiltonian' 

H = (k 2 + dlf + 30| (x) - e - 2q 2 (d 2 x + k 2 ) + q\ (22) 

The q = mode corresponds to an uniform displacement 
of the interface wherefore w 9= o is fixed to zero. Since 
■00 = this condition corresponds to 



= ect>' L - 30 



b' L -2^-^\ (23) 



which is satisfied because Eq.(|23|) is the derivative of 
Eq.pl). 
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FIG. 5. We here sketch (a) the unperturbed lamellar phase 
and (b) a modulated lamellar phase. The perturbation 
<j> — A(t) cos(qy)ifj q introduce spatial variations in the pat- 
tern at the length scale L — 2n/q. 



In the limit e — > we have <pL — * 0, so for e 
the eigenvalue equation reads —u> q ip q — — (fcg 







.,2\2 



2{ka — q 2 )ip q — ip' q ' ■ Remembering that ip q is periodic in x 
with wavevector kg, this implies uj q — q A . For small e we 



expect u> q 



aq 2 + q 4 7 



where a is a constant. We can ver- 



ify this assertion and determine the value of a by treat- 
ing the q-dependent part of H as a perturbation, and 
calculate the first-order correction to the ground state 
eigenvalue using standard perturbation techniques. The 
unperturbed eigenfunction for H is "0o = 4>' L with eigen- 
value zero. Obviously the g 4 -term in the perturbation 
just gives a contribution q 4 to the eigenvalue. The 0(q 2 ) 
contribution is 

Jdx^idl + k 2 )^ 



aq 2 = -2q 2 



Ux^y 



2q< 



JdxWl) 2 



J dx ( 



and with 0^ expanded as previously [Eq.(|17|)] we have 



e /256. Thus the relaxation rate is 



256 
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(24) 



We notice that the coefficient a here assumes the same 
value as determined above by the projection operator 
method, and that a dimensional analysis of Eq. (^4|) thus 
predicts the same crossover behaviour as did the analysis 
of Eq.(16) for the interface velocity. Furthermore, we 
notice that the formal expression for the coefficient a is 
identical to that obtained from the projection operator 
method, as may be seen from an integration by parts, 
i.e., the result holds generally, not just to the order given 
by the expansion Eq. (Jl7|) . 

Due to the more complicated structure of the DC- 
equation [Eq.(^)] a similar analysis of interfacial relax- 
ation in diblock copolymers has not yet proved possible. 
The main difficulty is that the 'Hamiltonian' operator H 
for this case is not self-adjoint, even for q = 0, with the 
result that a perturbative calculation of uj q requires not 
only the null eigenfunction 4>' L (x) of the q — operator 
Hq, but also the null eigenfunction of the adjoint opera- 
tor Hq, which we have so far been unable to determine. 



G 



V. SUMMARY AND DISCUSSION 

By numerical investigations we have found evidence 
of identical coarsening dynamics for the lamellar phase 
of the Swift-Hohenberg and diblock copolymer systems. 
This suggest that both systems belong to the same uni- 
versality class. We have extracted temperature depen- 
dent dynamical scaling exponents for the characteristic 
length scale partly by computing the ordinary structure 
factor and partly by computing a correlation function 
[Eq.(^|)] of the director field. Surprisingly the two meth- 
ods yield different scaling exponents indicating that the 
scaling phenomenon in question is non-trivial. We have 
no good understanding of the reasons for this discrep- 
ancy, but it should be noticed that the length scale ex- 
tracted from the structure factor does not have the same 
immediate geometrical interpretation as has the length 
scale extracted from the director-field correlation func- 
tion. Furthermore, the fact that the length scale, L nn , 
arising from the director field correlation function scales 
with the same growth exponents as the energy suggests 
that L nn is the physically important length scale in the 
system. 

Theoretically, by considering how curved interfaces re- 
lax we have demonstrated that the projection operator 
method, when applied to either of the two systems, re- 
sults in the same scaling exponents. This finding sup- 
ports the suggestion from our numerical work that the 
coarsening dynamics of the Swift-Hohenberg and di- 
block copolymer systems belong to the same universality 
class. However, the theoretical analysis applies only to 
defect free systems and does not explain the observed 
temperature-dependence of the growth exponents. A 
thorough understanding of the coarsening phenomenon 
here considered requires a theoretical treatment which 
successfully includes the simultaneous effects of both in- 
terfacial relaxation and defect-defect interactions. 
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